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Abstract. We present a theory for spin-polarized transport through a generic 
organic polymer connected to ferromagnetic leads with arbitrary angle 9 between 
their magnetization directions, taking into account the polaron and bipolaron states as 
effective charge and spin carriers. Within a diffusive description of polaron-bipolaron 
transport including polaron-bipolaron conversion, we find that the bipolaron density 
depends on the angle 9. This is remarkable, given the fact that bipolarons are spinless 
quasiparticles, and opens a new way to probe spin accumulation in organic polymers. 
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1. Introduction 

Recent years have witnessed significant advances in organic electronics, with interesting 
new fundamental insights and the prospect of new applications and devices functioning 
at room temperature [1,2]. A particularly interesting aspect comes from the spin 
degree of freedom, leading to "plastic spintronics" [3]. Organic materials such as 
polymers may be superior to inorganic semiconductor devices because of their small 
spin-orbit and hyperfine couplings, in principle allowing for very long spin coherence 
times. Moreover, the ease of fabrication and low-temperature processing of organic 
materials is very attractive for possible applications. Spin transport through tt- 
conjugated semiconducting organic polymers has consequently been studied in a number 
of recent experiments, and evidence for spin-polarized current injection and giant 
magnetoresistance in organic spin valves [4-7] as well as spin-dependent optical effects 
[8, 9] have been reported. 

Besides relevance for applications, the unconventional electronic properties of 
conducting polymers pose interesting fundamental questions. In undoped trans- 
polyacetylene, the charge and spin carriers are known to be soliton-like excitations, which 
are characterized by nontrivial spin-charge relations reficcting electron fractionalization 
[10]. This raises the possibility of unconventional spin-transport properties in undoped 
trans-polyacctylcnc. On the other hand, for basically all doped (nondcgcncratc) 
polymers, it has been established that the dominant charge and spin carriers at low 
energy scales well below the mean-field Peierls gap A correspond to polarons and 
hipolarons [1,10,11], whereas solitons can safely be ignored. As the polaron carries 
spin 1/2 like an ordinary electron and the bipolaron is spinless, spin current can only 
be carried by the polaron. Nevertheless, as we show below, the bipolaron density is 
affected by spin-polarized transport and can serve as a tool to detect the latter. 

In this work, we discuss spin transport through doped organic polymers, where 
polarons and bipolarons are the relevant charge carriers. In a typical two-terminal 
geometry (transport along the x axis), the organic polymer is contacted at ,t = and 
X = L hy two ferromagnetic (FM) metallic electrodes, where L is the length of the 
polymer. The left (right) electrode is characterized by a magnetization unit vector 

("^i?), with the angle 6 between them, • ttir — cos 9. We do not attempt 
a microscopic modelling of the interface between a FM electrode and the organic 
polymer, but follow the arguments of refs. [11-14], where it has been established that 
carriers injected into the polymer tunnel predominantly into polaron states close to 
the contact. We therefore impose the boundary condition that no bipolaron states 
near the boundaries (at a; = and x = L) arc filled by the injected current. Both 
contacts can then be completely described by spin- dependent conductances and 
G^, which take into account the spin-dependent density of states in the FM and the 
(disorder-averaged) matrix elements for tunneling into polaron states [15]. Moreover, 
for noncoUinear magnetizations (0 < ^ < tt), one also has to include the complex- valued 
mixing conductance reflecting boundary exchange processes [16-18]. 
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Transport in the polymer itself has so far been modelled either numerically, using 
lattice simulations of charge transport [19-22], or analytically, using simple master 
equations [23] or drift- diffusion models. The latter approaches have also been applied 
to spin transport [24-27]. Here we use the network theory of ref. [15,16] combined 
with a diffusive model to obtain spin-transport properties of a doped organic polymer 
sandwiched between two FM electrodes with noncoUinear magnetization directions 
(arbitrary 9). In the absence of bipolarons and for very high temperatures, this problem 
has been studied in ref. [26]. Here we present a generalization including the polaron- 
bipolaron conversion process, and also study the low-temperature quantum-degenerate 
limit. 



2. Model 

The energy-dependent polaron distribution function fp{x,e) at location < x < L 
can be decomposed into a spin- independent scalar part fo{x, e) and a spin-polarization 
vector f{x, e), 

fp{x,e) = foao + f ■ (T, (1) 

with the standard Pauli matrices cij in spin space; o"o is the unit matrix, and we assume 
homogeneity in the transverse direction. Note that a polaron has charge e and spin 
1/2. Another important charge carrier in organic polymers is the spinless bipolaron, 
with charge 2e and the scalar distribution function fBp{x,e) [1,10]. With the average 
density of states p(e), we introduce normalized densities by integrating the distribution 
functions over energy, 

hp{x) = J d€p{€)fp{x,€) = no{x) + n{x) ■ <T, 

nBp{x)^ J dep{e)fBp{x,e). (2) 

These densities are defined relative to an equilibrium reference value, and refiect 
nonequilibrium charge and spin accumulation in the polymer. Since our model does 
not include the quasiparticle states outside the mean-field gap A, but only retains the 
polaron and bipolaron states inside the gap, we choose p(e) ~ A~^©(A — |e|). 

In typical organic polymers, disorder is present and implies diffusive transport for 
both polarons and bipolarons, with the respective diffusion constants Dp and Dbp- The 
equations of motion for np{x,t) and nBp{x,t) are thus 

dtUp ^ Dpdlhp — r^f^ {hp — rioao) + i[h • cr,hp]- — Spao, (3) 
dtriBP = DBpdlriBP + Sp, (4) 

where Tgf is the polaron spin-relaxation time and Sp models conversion processes 
between polarons and bipolarons [27], 

Sp{x) = k (nl — n^) — bupp- (5) 
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The parameter k describes the local recombination rate for two polarons of opposite 
spin forming a bipolaron, while b comes from the reverse process, where a bipolaron 
decomposes into two polarons of opposite spin. The spin-precession term in (3) comes 
from an applied homogeneous magnetic field, where h = gfisB/h. We are interested 
in the steady-state case, where dthp = dtfiBP = in (3) and (4). According to Pick's 
law, the stationary spin- dependent particle current in the polymer is then encoded in 
the 2x2 matrix (in spin space) 

J{x) = -Dpd^hp{x) - DBpd^nBp{x)ao. (6) 

Equation (3) yields a decoupled equation for the spin polarization vector. 



Dpdln{x) = I r;/ -h, 

h"}i h/'j' 'T~^r 



■n{x). (7) 

/ 

Given the solution to (7), by taking the scalar part of (3) and combining it with (4), 
the bipolaron density is determined by 

nBp{x) = --^ (0^+V + no(x)) , (8) 

with two integration constants O and V. The only nontrivial equation that needs to be 
solved is given by 

.2^ ,.,2 , bDp fO 



Dpdino = k{ni - n') + [-x + V + nA. (9) 

J-' BP \^ J 

As discussed above, we impose the boundary condition 

nBp(O) = nBp{L) = 0, (10) 

since tunneling into the polymer involves only polaron states. With (8), this implies 
boundary conditions for (9), 

no(0) = -P, no(L) = -(P + 0). (11) 

In order to solve (7), we need six additional integration constants. Wc therefore have to 
specify boundary conditions reflecting spin and charge current continuity at the contacts 
to the left and right FMs. The FMs are taken as reservoirs with identical temperature 
T and chemical potentials iiliri where the applied voltage is eV — iiL — i^R- As before, 
we introduce (normalized) densities, 

y dep{e)nFie - (Il/r), (12) 

with the Fermi function npie) = l/[e'^/^^'^ + 1]. Boundary conditions then follow by 
relating the current (6) aX, x — Q {x — L) io the injected current at the left (right) 
interface [16], 

J(0) = - ^ G'^ul{hp{Q)-nl^a^)ul - (^^^^^^^(o)^ + h.c.) , (13) 
J(L) = G^uUML) - nr<^o)^S + {g^^uIML)u\, + h.c.) . (14) 



^L/R 
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Note that (10) implies that bipolarons do not enter this boundary condition. The 
matrices "uj ^ = ^{^ + o'7fiL,R-cr) project the spin direction a =|, [— +, — in the polymer 
onto the respective FM magnetization direction. For simplicity, we assumed identical 
spin-polarized {G\G^) and mixing (G^^) conductances for both contacts. They must 
obey ReG^^ > {G^ + G^)/2 [16]. The 2 x 2 matrix equations (13) and (14) allow to 
determine the eight integration constants, and thereby yield the spin-polarized current 
through the system for arbitrary 9. Moreover, this gives access to the bipolaron density 
from (8) after solving (9). We stress that none of the eight integration constants depends 
on the parameters k and b in (5). 

From (6) and (8), we can immediately see that charge current Jc — DpO/L is 
conserved, 

J{x) ^ JcOo^ J s{x) ■ a, (15) 

and the spin current, Js{x) — —Dpdxn{x), follows from the solution of (7). Remarkably, 
both Jc and Js{x) are independent of the polaron- bipolaron transition rates k and b in 
(5), and the spin-dependent current alone cannot detect the presence of bipolarons in 

the polymer. Nevertheless, as we show below, the bipolaron density nBp{x), which is 
induced by the nonequilibrium spin accumulation in the polymer, is sensitive to these 
rates. As a useful measure, we will employ the integrated density, 

N{9) ^ [ dx nBp{x]d). (16) 
The 6'-dependence of the bipolaron density is then encoded in the dimensionless quantity 

^^^^ - N{0)-N{ny ^^^^ 

By definition, this quantity interpolates between R{0) = and R{7t) = 1 as ^ is varied 
from the parallel to the antiparallel configuration. 

3. Collinear case: a readily solvable limit 

We first discuss a simple yet important limit, where a direct analytical solution can 
be obtained. This limit is defined by collinear magnetizations, mR — prhi with 
p = ± (parallel or antiparallel configuration) and = e^. Moreover, we consider the 
length of the polymer as short compared to the spin coherence length, L <^ ^yDpTsf, 
and put h = (no magnetic field). In that case, (7) has the general solution 
n{x) — —{Fx/L + G), with constant vectors F and G. For rhi, — — -^rhp, the 
boundary conditions (13) and (14) imply that the x and y components of both vectors 
vanish, and the spin current is conserved, 

n{x) = -eA'^x + g\ , Js = ^^e^. (18) 
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The four remaining integration constants {O, V, T ^ Q) readily follow by solving the 
boundary conditions (13) and (14) [16]. For the parallel {p — +) configuration, they are 

= (GT + 2Gp)(Gi + 2Gp)'^' ^'^^ 
while for the antiparallel case, we find J-'- = and 

0_ - 2(7^_ + Q^Qi^2Gp{G^ + GI) ' 
C = - (g^-G^)Gp eV 

^- GTCi + 2G'p(GT + Gi) 2 ' ^ ^ 

where /2 = (/^l + f^R)/'^ is the mean chemical potential and Gp = Dp/L. The charge 
current for the respective configuration is then Jc = GpO±, while the spin current is 
Js — GpT±ez- 

The remaining task is to solve (for given p = ±) the nonlinear equation (9) for 

no{x) under the boundary condition (11), using (18), (19) and (20). Since the transition 
rates k and h are known to be small [27], we use a perturbative iteration scheme and 
write 

no(x) = -^x-P + no(x). (21) 

For k — h — this Ansatz solves (9) under the correct boundary conditions when 
putting nQ{x) — 0; note that the bipolaron density is directly proportional to fio(a:), see 
(8). For small but finite /c, 6, we then insert (21) into (9) and linearize in fiQ. This yields 
a second-order differential equation for no(x), which needs to be solved under Dirichlet 
boundary conditions at a; = and x = L. The solution gives the bipolaron density for 
the parallel and antiparallel configuration in closed form, 



Dbp V ^ ^ 2 L 



+ T7T K (22) 



L2 12^ 

where the integration constant C± follows from the condition nBp{L) ~ 0. The 
integrated bipolaron density (16) is then given by 

N± = (vi -gl + - T^g^ + ^(^1 - ^l)) ■ (23) 

Note that ^ follows immediately from (19) and (20), indicating that the 
bipolaron density indeed is sensitive to the spin accumulation in the polymer. The 
bipolaron density (22) is shown in figure 1, taking parameters for sexithienyl as 
organic spacer [26]. One clearly observes a difference between the parallel and the 
antiparallel configuration. Although the current is not sensitive to the polaron-bipolaron 
transition rates k and 6, the bipolaron density is influenced by the nonequilibrium spin 
accumulation in the polymer. 
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x/L 



Figure 1. Bipolaron density nBp{x,6) for coUincar magnetizations, i.e., 6 = (solid 
curve) and 6 = n (dashed curve), as obtained from (22). We use a representative 
parameter set for hole transport [26]. In units where Gp = Dp/L = 1, the parameters 
are Dbp = 2L/3,GT = lO^,^^ = 10-2,fc = Dp/10, b = Dbp/10, A = 3.5, p- = -2, 
T = 0, and eV = 1. The inset shows y{x) = [n^"^(x) - n^^(x)]/n^"^(L/2) for the 
same curves. 
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Figure 2. Spin-accumulation sensitivity R{9), see (17), of the bipolaron density as a 
function of the magnetization tilt angle. Parameters are as in figure 1, additionally we 
set ReGTi = ImG^i = 5.1 x lO^. 

4. Noncollinear magnetization 

In the general case of arbitrary angle 9 between and rhji, one can solve the problem in 
an analogous manner but the equations become less transparent. The main difference 
is that now the mixing conductance G'^^ has to be taken into account. However, as 
reported previously [26], we find that the results are practically independent of the 
precise choice for G'^^. We find a smooth crossover between the limiting values for = 
and 9 = 7:, see (22), illustrated for the integrated bipolaron density (16) in figure 2. 



Spin transport and bipolaron density in organic polymers 



8 



5. Conclusions 

In this work, we have discussed spin transport in doped organic polymers, employing a 
diffusive description of polaron and bipolaron transport. In a two-terminal setup, where 
the polymer is sandwiched by (generally noncoUinear) ferromagnetic electrodes, the 
problem can be solved analytically by exploiting the smallness of the polaron-bipolaron 
transition rates k and b. While the spin-dependent current through the device turns out 
to be independent of k and b, the nonequilibrium bipolaron density is a sensitive probe 
of spin accumulation. The possibility to measure this density in optical-absorption 
experiments [28-30], e.g., by adapting charge- modulation techniques [31] to the two- 
terminal transport geometry considered here, may offer a novel way to probe spin 
accumulation in organic polymers. Such an optical method would be complementary to 
the usual magnetoresistance measurement of spin accumulation and could thus serve as 
another means to independently verify spin-injection efficiencies in organic polymers [32]. 

Our work generalizes previous studies where bipolarons were neglected [26] or 
only a single ferromagnet-polymer interface was considered [27]. We also treat 
the nonequilibrium situation due to an applied voltage self-consistently instead of 
postulating the existence of a uniform electric field [27]. Wc mention in passing that 
results from a recent Monte Carlo simulation [33] have elucidated the importance 
of bipolaronic effects for a nontraditional type of magnetoresistance that occurs in 
conducting polymers in the absence of magnetic contacts. Another recent theoretical 
study [34] on magnetoresistance in polymers with polaron and bipolaron carriers used a 
diffusive approach and magnetic contacts (FM-polymer-FM configuration). In contrast 
to our work, however, ref. [34] does not take into account conversion processes between 
polaron and bipolaron states, but simply assumes a constant density of bipolarons 
and includes this into the transport calculations. Surprisingly, a dependence of the 
magnetoresistance on the ratio of bipolarons and polarons is reported [34], whereas we 
find the spin-polarized current to be independent of the bipolaron formation rate. Our 
finding can be traced back to the well-established [13, 14] suppression of tunneling into 
bipolaron states near the interface with a FM electrode. This feature is ignored when 
simply asssuming a constant bipolaron density. 
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